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Two new luminosity functions of galaxies can be built starting from 
three and four parameter generalized gamma distributions. In the astro- 
physical conversion, the number of parameters increases by one, due to the 
addition of the overall density of galaxies. A third new galaxy luminosity 
function is built starting from a three parameter generalized gamma distri- 
bution for the mass of galaxies once a simple nonlinear relationship between 
mass and luminosity is assumed; in this case the number of parameters is 
five because the overall density of galaxies and a parameter that regulates 
mass and luminosity are added. The three new galaxy luminosity functions 
were tested on the Sloan Digital Sky Survey (SDSS) in five different bands; 
the results always produce a "better fit" than the Schechter function. The 
formalism that has been developed allows to analyze the Schechter function 
with a transformation of location. A test between theoretical and observed 
number of galaxies as a function of redshift was done on data extracted 
from a two-degree field galaxy redshift survey. 

PACS numbers: 02.50.Cw Probability theory; 98.62.Ve Statistical and correl- 
ative studies of properties (luminosity and mass functions; mass-to-light ratio ) 



1. Introduction 

The luminosity function of galaxies &(L)dL (LF) is the number of galax- 
ies per unit volume, Mpc s , whose luminosity is comprised between L and 
L + dl, see Section 3.7 in pQ or Section 1.12 in [2]. The luminosity, L, has 
physical units of Watt Hz -1 and therefore the band being considered should 
always be specified. In our case, we selected the data of the Sloan Digital 
Sky Survey (SDSS) which has five bands u* (A = 3550A), g* (A = 4770A), 
r* (A = 6230A), i* (A = 7620A) and z* (A = 9130A) with A denoting the 
wavelength of the CCD camera, see [3]. 
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The Schechter LF of galaxies, see [3], is the more widely used function 

mdL = (^)(j ? reM-§-)dL , (1) 

where a sets the slope for low values of L, L* is the characteristic luminosity 
and $* is the normalization. In the formula above, L* characterizes the 
break of the LF. As an example, L > 0.1L* defines a giant galaxy , see the 
text after formula (1.18) in [5]. An astronomical form of equation (pQ) can 
be deduced by introducing the distribution in absolute magnitude 

<$>{M)dM = (0.4ZnlO)$*10 a4(a+1)(M *- M) 

x exp(-10 a4 ( M *- M VM , (2) 

where M* is the characteristic magnitude as derived from the data. 

This function has three parameters that can be found by fitting the data. 
Over the years, many modifications have been made to the Schechter LF in 
order to improve its fit: we report two of them. When the fit of the rich 
clusters LF is not satisfactory a two-component Schechter-like function is 
introduced, see [6] 



<£* L L 
L max > L > L Dwarf : ${L)dL = (jj)(j-;T exp(-— )dL 



• & Dwarf / L 



(3) 



Lowarf >L> L mm : $(L)dL = (-^L)(_ ^ 

where 



L* Lz) war f 



fL Dwarf \a / Lj^warf \ 
Dwarf = $ ( L , ) exp( j^L) . 

This two-component function defined between the maximum luminosity, 
Lmax, and the minimum luminosity, L m j n , has five parameters because two 
additional parameters have been added: Lowarf which represents the mag- 
nitude where dwarfs first dominate over giants and otDwarf which regulates 
the faint slope parameter for the dwarf population. 

Another LF introduced in order to fit the case of extremely low luminos- 
ity galaxies is the double Schechter function with five parameters, see [TJ: 



${L)dL = ^exp(-L/L* 



L\ ai /L^" 2 



(4) 



where the parameters and a which characterize the Schechter function 
have been doubled in <j>^\ and 0*^. The strong dependence of LF on dif- 
ferent environments such as voids, superclusters and supercluster cores was 
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analyzed by [8] with 

F(L)dL oc (L/L*) a (l + (L/Ly^'^^^L/L*), (5) 

where a is the exponent at low luminosities (L/L*) <C 1, 5 is the exponent 
at high luminosities (L/L*) S> 1, 7 is a parameter of transition between the 
two power laws, and L* is the characteristic luminosity. The previous LFs 
leave a series of questions unanswered or partially answered: 

• What is the function of introducing a transformation of location in 
the LF? 

• Is it possible to model the data with just a single LF? 

• Is it possible to deduce a LF for galaxies starting from the mass dis- 
tribution of galaxies? 

• Is it possible to improve the Schechter function by introducing a trans- 
formation of location? 

• Does the new LF match the observed behavior of the number of galax- 
ies for a given solid angle and flux as a function of the redshift? 

In Section [21 this paper explores how a generalized gamma distribution can 
model two galaxy LFs. The method which allows the LF to be deduced for 
galaxies starting from a mass distribution as given by a generalized gamma is 
presented in Section[3j Section|3]reports the analytical and numerical results 
on the Schechter function with transformation of location. In Section the 
redshift dependence of the Schechter function and one of the four new LFs 
are compared with data from the two-degree Field Galaxy Redshift Survey 
in the 2dFGRS catalogue. 

2. The generalized gamma distribution 

The generalized gamma distribution can be represented with three pa- 
rameters, see [9], [10] or four parameters, see We will explore both cases 
in the following. In order to make a comparison between our LF and the 
Schechter LF, we first down-loaded the data of the LF of galaxies in the five 
bands of SDSS adopted in [12] : they are available at: http:/ /cosmo. nyu.edu/blanton/lf.html 
In the previous paper, |12j . the basic assumption was to consider a Friedmann- 
Robertson- Walker cosmological world model with matter density fio = 0.3, 
vacuum pressure Q\ = 0.7 and Hubble constant Hq = 100 h km s _1 Mpc -1 
with h = 1. The data contain the absolute magnitude, the value of the LF 
for that magnitude and the error of the LF. 
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The LF of galaxies as obtained from astronomical observations ranges 
in magnitude from a minimum value, M m i n , to a maximum value, M max ; 
details can be found in [13] and [13] • A nonlinear fit through the Levenberg- 
Marquardt method (subroutine MRQMIN in [15]) allows the determination 
of the parameters, but the first derivative of the LF with respect to the 
unknown parameters should be provided. The merit function x 2 can be 
computed as 

j — l ^LFastr 

where n is number of data and the two indices theo and astr stand for 
theoretical and astronomical, respectively. 

Particular attention should be paid to the number of unknown parame- 
ters in the LF: three for the Schechter function (formula ([2])), four for for- 
mula (|2ip . four for the Schechter function with transformation of location 
(formula (|30|) ). five as represented by formula (|28|) and five for formula (|16|) . 
A reduced merit function x 2 e d can be computed as 

Xred = X 2 /NF , (7) 

where NF = n — k, n is the number of data and k is the number of param- 
eters. The Akaike information criterion (AIC), see |16] , is defined as 

AIC = 2k- 2ln(L) , (8) 

where L is the likelihood function and k the number of free parameters in the 
model. We assume a Gaussian distribution for the errors and the likelihood 
function can be derived from the x 2 statistic L cx exp(— ^-) where x 2 has 
been computed by equation ([6]), see [17], [18]. Now AIC becomes 

AIC = 2k + X 2 ■ (9) 

The Bayesian information criterion (BIC), see |19j . is 

BIC = k ln{n) - 2ln(L) , (10) 

where L is the likelihood function, k the number of free parameters in the 
model and n the number of observations. The phrase "better fit" used in 
the following means that the three statistical indicators : x 2 ; AIC and BIC 
are smaller for the considered LF than for the Schechter function . 
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2.1. The generalized gamma distribution with five parameters 

The starting point is the probability density function (in the following 
PDF) named generalized gamma that we report exactly as in 

G(x;a,b,c,k) = K 6 ; , 11 

br (c) 

where T is the gamma function, a is the location parameter, b is the scale 
parameter, c and k are two shape parameters. The number of parameters 
is four and the astrophysical version of the previous PDF can be obtained 
by inserting a = L a , x = L and b = L*: 



*(L;L*,c,k,**) = ** K L * ^ , (12) 

where is a normalization factor which defines the overall density of galax- 
ies, expressed as a number per cubic Mpc. The mathematical range of exis- 
tence is L a < L < oo and the number of parameters is five because \P* has 
been added. The averaged luminosity, (L), is: 



L*T(±&) +L a T(c) 



< L > = Tic) 1 : ' ' < U > 



L=( C ^-Y L* + L a . (14) 



and the mode is at 

T- = i 

V k 

The relationships connecting the absolute magnitude M, M* and M a of a 
galaxy to its luminosity are: 

10 0.4(M« Q -M) -^* = 1 q0.4(M, q -M*) = j_QO.4(Af» -Af„) 



L*0 LjfQ LjfQ 

(15) 

where M*0 is the absolute magnitude of the sun in the considered band. 
As an example, the SDSS bands have M* Q = 4.48 in z* and M*0= 6.32 in 
u* , see [20]. A more convenient form of the LF in terms of the absolute 
magnitude M is: 



0.4 A; In (10) f i °- 4M *-°- 4M - 10 
■%>(M)dM = <P* * 



rb— 1 

0AM*-0AM a x 



r( c ) 

( 10 0.4M*-0.4M_ 10 0.4M*-0.4M o ) fc .4M*-0.4M 



10 U - 4M ~ UAM dM .(16) 



G 
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Table 1. Parameters of fits to LF in SDSS Galaxies of five parameter function 
represented by formula (1161) . 



parameter 


u 


9 


r 


i 


z 


M a 


-14.25 


-15.55 


-15.77 


-16.52 


-17.39 


M* - 5 log 10 h 


-17.22 


-18.29 


-20 


-20.01 


-20.25 


^*[h 3 Mpc- 3 ] 


0.078 


0.077 


0.11 


0.05 


0.055 


c 


0.61 


0.30 


0.19 


0.47 


0.49 


k 


0.82 


0.92 


0.87 


0.78 


0.76 


N 


483 


599 


674 


709 


740 




-20.65 


-22.09 


-22.94 


-23.42 


-23.73 




-15.78 


-16.32 


-16.30 


-17.21 


-17.48 


K 


282 


736 


1766 


1726 


2136 


X-red 


0.59 


1.24 


2.64 


2.45 


2.90 


AIC k = 5 


292 


746 


1776 


1736 


2146 


BIC k = 5 


313 


768 


1799 


1759 


2169 


X 2 — Schechter 


330 


753 


2260 


2282 


3245 


Xred ~ Schechter 


0.689 


1.263 


3.368 


3.232 


4.403 


X 2 - Blanton 2003 


341 


756 


2276 


2283 


3262 


AIC k = 3 — Schechter 


336 


759 


2266 


2288 


3253 


BIC k = 3 — Schechter 


349 


772 


2279 


2302 


3265 



The mode when expressed in magnitude is at 

M = -1.085 In ( (^^) 1 e" - 921 M * + e" ' 921 M ^ . (17) 

This data-oriented function contains the five parameters c, k, M*, M a 
and which can be derived from the operation of fitting observational 
data. The results are reported in Table Q] together with the number of 
elements N belonging to the sample, M m i n and M max of the sample, the 
merit function x 2 an d Xred' * ne X 2 an d X 2 - e d °^ ^ ne Schechter function as 
computed by us and the x 2 of the Schechter function as computed by [P2"] . 

The Schechter function, the new five parameter function represented by 
formula (fT6j) and the data are reported in Figure [lj Figure [21 Figure [3j 
Figure [U Figure El Figure El Figure [71 Figure El Figure El and Figure flQl 
where bands u* , g*, r*, i* and z* are considered. 

The flexiblity of the generalized gamma with five parameters to fit a 
sudden increase in the low luminosity region is clearly visible in Figure [TUl 
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-16 
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Fig. 1. The luminosity function data of SDSS(u*) arc represented with error bars. 
The continuous line fit represents our LF (116[) and the dotted line represents the 
Schechter function. 

2.2. The generalized gamma distribution with four parameters 
We can start from equation (|lip . inserting a = 0, x = L and b = L*: 



The mathematical range of existence is < L < oo and the number of 
parameters is four because a = and \P* have been added. The averaged 
luminosity is 




(18) 




(19) 



and the mode is at 




(20) 



The magnitude version of the LF is 



**0.4 In (10) fclO-°- 4cfe ( M - M *)e- 10 



0.4 (M- M*)k 



^(M)dM 



dM 



(21) 
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absolute magnitude M 



Fig. 2. The luminosity function data of SDSS(u*) in the low luminosity region are 
represented by the error bar. The continuous line fit represents our LF (|16[) . the 
dotted line represents the Schechter function as given by our data and the dashed 
line represents the Schechter function as given by |12j . 

The mode when expressed in magnitude is at 



This function contains the four parameters c, kk, M* and \&* and the 
numerical values obtained are reported in Table [2j 



We assume that the masses of galaxies, Ai, are distributed as a gener- 
alized gamma. We can start from equation (jlip inserting a = , x = Ai 




(22) 



3. The generalized gamma distribution for mass 



and b = M* 



y{M;M*,c,k : ^*) = ^ 



(23) 



A4*r ( c ) 



This is a generalized gamma distribution with a scale parameter M. * , c and k 
are shape parameters . The mathematical range of existence is < Ai < oo 
and the number of parameters is four because ^* has been added. The 
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Table 2. Parameters for fits to LF in SDSS Galaxies with the four parameter 
function represented by formula (|21[) . 



parameter 


u* 


9* 


r* 


i* 


z* 


M* - 5 log 10 h 


-17.34 


-19.43 


-20.25 


-20.29 


-20.76 


$>*[h 3 Mpc~ 3 ] 


0.101 


0.231 


0.65 


0.097 


0.2 


c 


0.47 


0.08 


0.02 


0.247 


0.10 


k 


0.842 


1.019 


0.93 


0.839 


0.864 




283 


747 


2188 


1866 


2915 




0.591 


1.256 


3.266 


2.648 


3.991 


AIC k = 4 


291 


755 


2196 


1874 


2921 


BIC k = A 


307 


773 


2214 


1893 


2921 


X 2 Schechter 


330 


753 


2260 


2282 


3245 


Xred ~ Schechter 


0.689 


1.263 


3.368 


3.232 


4.403 


X 2 Blanton 2003 


341 


756 


2276 


2283 


3262 



average value is 



r 



l±ck\ 

(M) = V*M* l/V • (24) 

r ( c ) 



10 
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Fig. 4. The luminosity function data of SDSS(g*) in the low luminosity region are 
represented by the error bar. The continuous line fit represents our LF (|16[) . the 
dotted line represents the Schechter function as given by our data and the dashed 
line represents the Schechter function as given by |12j . 

We now follow the pattern fixed in [21J where the PDF representing the 
mass of galaxies was assumed to be a Kiang function, see [22]. The mass- 
luminosity relationship is assumed to follow a power law having the form 



where 1/d is an exponent which connects the mass to luminosity. The PDF 
([23l ) is therefore transformed into the following form: 



The mathematical range of existence is < L < oo and the averaged 
luminosity, (L), is 




(25) 



^(L)dL = * 



cfc— 1 k 1 



(26) 



dF{c)L 




(27) 
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absolute magnitude M 



Fig. 5. The luminosity function data of SDSS(r*) are represented with error bars. 
The continuous line fit represents our LF (116[) and the dotted line represents the 
Schechter function. 



The magnitude version of the previous LF is: 



^{M)dM = V* 



0.4 ak 10~ - 4 



cak (- M*+M) _ 1q -0.4 



d e 



(- M*+M)ak 



In (10) 



T(c)d 



(28) 



The mode expressed in magnitude is at 



M = M* - 1.0857 



In (c) d 



(29) 



This function contains the five parameters c, k, d, M* and and the 
numerical values obtained are reported in Table O 

From a careful analysis of TableO it is possible to conclude that 1.015 < 
d < 1.024. Because we have assumed that L oc Ai d , the nonlinear relation- 
ship between mass and luminosity is a weak effect. A weak dependence of 
the parameter that regulates mass and luminosity has also been found in 
|23| when the PDF which represents the mass of galaxies is the product of 
two gamma variates with argument 2. 
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-16.6 



absolute magnitude M 



Fig. 6. The luminosity function data of SDSS(r*) in the low luminosity region are 
represented with error bars. The continuous line ht represents our LF (|16|) . the 
dotted line represents the Schechter function as given by our data and the dashed 
line represents the Schechter function as given by |12) . 



Table 3. Parameters of fits to LF in SDSS Galaxies of the five parameter function 
represented by formula (|28p deduced from the mass-luminosity relationship. 



parameter 


u* 


9* 


r* 


i* 


z* 


M* - 5 log 10 h 


-17.34 


-19.39 


-20.18 


-20.43 


-20.53 


f*[/i 3 Mpc~ 3 ] 


0.101 


0.19 


0.33 


0.12 


0.11 


c 


0.47 


0.1 


0.05 


0.17 


0.22 


k 


0.855 


1.029 


0.93 


0.88 


0.824 


d 


1.015 


1.024 


1.023 


1.018 


1.017 


x\ 


283 


751 


2217 


1896 


3009 


X-red 


0.592 


1.265 


3.313 


2.693 


3.313 


AIC k = 3 


289 


757 


2223 


1902 


3015 


BIC k = 3 


314 


783 


2249 


1929 


3042 


X 2 Schechter 


330 


753 


2260 


2282 


3245 


Xred ~ Schechter 


0.689 


1.263 


3.368 


3.232 


4.403 


X 2 Blanton 2003 


341 


756 


2276 


2283 


3262 
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I sample 



...■'■"''"J-'k 




-23 



-22 -21 -20 _ -19 

absolute magnitude M 



-18 



Fig. 7. The luminosity function data of SDSS(i*) are represented with error bars. 
The continuous line fit represents our LF (I16[) and the dotted line represents the 
Schechter function. 



4. The Schechter function with transformation of location 

The Schechter LF of galaxies once the location L a is introduced is: 

HL)dL = (^)( i£ T ^)° eM JJ ^)dL , (30) 

and is characterized by four parameters. 
The average luminosity, (L), is 

(L) = <P*(L*T(2 + a)+L a T(a + l)) , 

and the mode is at 

L = a L* + L a , 



or 



M = -1.085 In (a e" a921 M * + e" a921 M " 



The magnitude version of this LF is 



<f>(M)dM = 0.4$* (10 



MAM*-0AM 



10 



. 1Q 0.4 M* -0.4 M +10 0.4 M* -0.4 M a QA M » _ Q 4 M 



0.4 A/*— 0.4 M 



In (10) 



(31) 

(32) 
(33) 



(34) 
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-17.6 



17.4 



-17.2 



absolute magnitude M 



Fig. 8. The luminosity function data of SDSS(i*) in the low luminosity region are 
represented with error bars. The continuous line fit represents our LF (jT6j) , the 
dotted line represents the Schechter function as given by our data and the dashed 
line represents the Schechter function as given by |12j . 



Table 4. Parameters of fits to LF in SDSS Galaxies of the four parameter Schechter 
function with location as represented by formula (1341) . 



parameter 


u* 


9* 


r* 


i* 


z* 


M* - 5 log 10 h 


-17.92 


-19.37 


-20.41 


-20.83 


-21.15 


$*[h 3 Mpc~ 3 ] 


0.03 


0.021 


0.015 


0.014 


0.013 


a 


-0.9 


-0.84 


-0.99 


-1.02 


-1.01 


M a [mags] 


-12.1 


-14.87 


-15.14 


-12.72 


-16.39 


K 


333 


723 


2108 


2453 


3133 




0.69 


1.21 


3.14 


3.47 


4.25 


AIC k = 4 


341 


731 


2116 


2451 


3141 


BIC k = 4 


358 


748 


2134 


2479 


3159 


X 2 Schechter 


330 


753 


2260 


2282 


3245 


Xred ~ Schechter 


0.689 


1.263 


3.368 


3.232 


4.403 


X 2 Blanton 2003 


341 


756 


2276 


2283 


3262 



This function contains the four parameters a, M a , M* and and the 
numerical values obtained are reported in Table HI 
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Fig. 9. The luminosity function data of SDSS(z*) are represented with error bars. 
The continuous line fit represents our LF (116[) and the dotted line represents the 
Schechter function. 

From a careful analysis of Table 01 it is possible to conclude that in three 
cases out of five analyzed, the Schechter LF with transformation of location 
produces a "better fit" than the standard version. 



Some useful formulae connected with the Schechter LF, see formula ([I]), 
in a Euclidean, non-relativistic and homogeneous universe are reviewed; by 
analogy, new formulae for the generalized gamma distribution with four 
parameters, see formula (fTHj) . are derived. 



The radiation flux, /, is introduced 



where r represents the distance of the galaxy. The joint distribution in 
redshift, z, and / adopting the Schechter LF, see formula (1.104) in [23] or 
formula (1.117) in [2J, is: 



5. The radial distribution of galaxies 



5.1. Dependence of the LF on z 



dVtdzdf 



dN 




(36) 
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absolute magnitude 



Fig. 10. The luminosity function data of SDSS(z*) in the low luminosity region are 
represented with error bars. The continuous line fit represents our LF (|16|) . the 
dotted line represents the Schechter function as given by our data and the dashed 
line represents the Schechter function as given by |12j . 



where dQ, dz and df represent the differential of the solid angle, the redshift 
and the flux, respectively. The critical value of z, z cr %u is 

2 _ ( n 7) 
mt ~ 47T/4 ' 1 ' 

where cl is the velocity of light and Hq is the Hubble constant. 

The number of galaxies, Ng(z, f m i n , fmax) comprised between a mini- 
mum value of flux, f m i n , and maximum value of flux f max , for the Schechter 
LF can be computed through the following integral: 

C fmax nj r- . 

Ns(z) = / 4vr(-f . (38) 

J fmin 11 z crit 

This integral does not have an analytical solution and we must perform a 
numerical integration. 

The number of galaxies for the Schechter LF in z and / as given by 
formula (|36p has a maximum at 



Zpos—max — %critV & H~ 2 , 



(39) 
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which can be re-expressed as 

_ ^2+WlO°- 4M *e- a4M *ff 

Zpos—max — n , — i-p > x™) 

where M*q is the reference magnitude of the sun in the bandpass under 
consideration. For the sake of clarity, we report the generalized gamma LF 
of galaxies with four parameters (equation (JTHJ) ) , in the following LF4 



.*(*) 



ck— 1 _i l \ k 
e V l* ) 



*(L;L*,c,k,**) = M>* v 'j—- . (41) 

The joint distribution in z and / of the generalized gamma LF4 of galaxies 
is 

The number of galaxies, Nlfa{z-, fmini fmax) of the LF4 comprised between 
f m in and fmax-, can be computed through the following integral: 

N LF4 (z)= f maX A*{^-fz^{^-)df , (43) 

J fmin -"0 Z crit 

and in this case again a numerical integration must be performed. 

The number of galaxies of the LF4 as given by formula (|42p has a max- 
imum at Zpos-max where 

jy 2 ln(l + cfc)-ln(fc) 

Zpos—max — 6 * Z cr ii , (,^4J 

which can be re-expressed as 

_ e l/2 ln(1+Ck r n(k) V 10 OAM* e -OAM* Hok 
Zpos—max — „ . — /-tt • V^^V 

The formulae previously derived are now tested on the 2dFGRS catalogue 
available at the web site: |http: / / msowww.anu. edu.au/2dFGRS7| The 2dF- 
GRS catalogue contains redshifts for 221,414 galaxies brighter than a mag- 
nitude limit of bJmag=19.45. The galaxies cover an area of approximately 
1500 square degrees and more details can be found in [25]. In particular 
we added together the file parent.ngp.txt which contains 145,652 entries for 
NGP strip sources and the file parent.sgp.txt which contains 204,490 entries 
for SGP strip sources. Once the heliocentric redshift had been selected, we 
processed 219,107 galaxies with 0.001 < z < 0.3. The parameters of the 



18 



MS PRINTED ON APRIL 28, 2010 

Table 5. Parameters of the Schechter function for 2dFGRS. 



parameter 2dFGRS 

M* - 5 log 10 h [mags] (-19.79 ± 0.04) + hLogh 

a -1.19 ±0.01 

h 3 [Mpc- 3 ] ((1.59 ± 0.1)10~ 2 ) x h 3 



Table 6. Parameters of LF4, formula ([32)1. 

based on the 2dFGRS data (triplets generated by the author) 

2dFGRS 

c 0.016 ±0.01 

M* — 51og 10 h[mags] -19.15 ±0.029 
**[/t 3 Mpc~ 3 ] 2.24 ±1.42 

k 0.82 ±0.06 



Schechter function for 2dFGRS are reported in Table [5) see first line of 
Table 3 in [26]. 

It is interesting to point out that other values for h different from 1 shift 
all the absolute magnitudes by 5 log 10 ^ an d change the number densities 
by a factor h 3 . 

Figure [TT] and Figure [12] report the number of observed galaxies in the 
2dFGRS catalogue for two different apparent magnitudes and the theoretical 
curves as represented by formula (|36p and formula (|42p . 

Due to the importance of the maximum as a function of z in the num- 
ber of galaxies, Figure [T3] reports the observed histograms in the 2dFGRS 
database and the theoretical curves as a function of magnitude. 

The total number of galaxies in the 2dFGRS database is reported in 
Figure [T3] as well as the theoretical curves as represented by the numerical 
integration of formula f|36j) and formula (|42p. 

Particular attention should be paid to the concept of limiting magnitude 
and to the corresponding completeness in absolute magnitude of the consid- 
ered catalogue as a function of redshift. The observable absolute magnitude 
as a function of the limiting apparent magnitude, mi, is 

M L = m L - 5 log 10 - 25 . (46) 
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Fig. 11. The galaxies of the 2dFGRS database with 16.77 < bJmag < 18.40 or 
4677 ^*^ < / < 21087 ^*g 2 (with bJmag representing the relative magnitude used 
in object selection), are isolated in order to represent a chosen value of m and then 
organized as frequency versus heliocentric redshift, (empty circles); the error bar 
is given by the square root of the frequency. The maximum in the frequencies of 
observed galaxies is at z = 0.085 when M*q = 5.33 and h=l . The theoretical curve 
generated by the Schechter function of luminosity (formula (|36l) and parameters as 
in column 2dFGRS of Table [S]) is drawn (full line). The theoretical curve generated 
by LF4, formula (|4"2"|) , and parameters as in column 2dFGRS of Tabled is drawn 
(dashed line); 8078 for the Schechter function and x 2 = 6654 for LF4. 

The interval covered by the LF of galaxies, AM, is defined as 

AM = M max - M min , (47) 

where M max and M m i n are the maximum and minimum absolute magni- 
tudes of the LF for the considered catalogue. The real observable interval 
in absolute magnitude, AMl, is 

AM L = M L - M min . (48) 

We can therefore introduce the range of observable absolute maximum mag- 
nitude expressed in percent, e(z), as 



(49) 
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Fig. 12. Galaxies in the 2dFGRS with 14.79 < bJmag < 15.02 or 105142-j^ < 
/ < 129757 ^*^2 ■ The maximum in the frequencies of observed galaxies is at 
z = 0.02, M* Q = 5.33, x 2 = 348 for the Schechter function (full line) and x 2 = 243 
for LF4 (dashed line). 

This is a number that represents the completeness of the sample and can be 
considered to be a version in percentage terms of the Malmquist bias for the 
apparent magnitude of a limited number of samples, see [271 EE]- Figure [TBI 
shows the behavior of the range in absolute magnitude in percent as a 
function of z when the limiting magnitude of 2dFGRS is m^=19.61 . From 
the previous figure, it is possible to conclude that the 2dFGRS catalogue is 
complete for z < 0.0442 . 

In this section we have adopted the absolute magnitude of the sun in 
the bj filter M* Q = 5.33, see [H1E2]. 

6. Summary 

We have presented three new LFs of galaxies based on two different ver- 
sions of the generalized gamma distribution. Before continuing, we should 
spend some time discussing the high values of x 2 obtained in Tables [TJ [2j [3] 
and[H We quote the following phrase from web site http: / /cosmo. nyu.edu/blanton/lf.html: 
"Thus, when comparing to theory it would be totally inappropriate to rely 
on x 2 from these errors. Such a x 2 will always be very high even for very 
successful theories." From the previous statement we can deduce that the 
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Fig. 13. Value of z at which the number of galaxies in the 2dFGRS database is 
maximum as a function of the apparent magnitude bJmag (stars) , theoretical curve 
of the maximum for the Schechter function as represented by formula (|40p (full 
line) and theoretical curve of the maximum for LF4 as represented by formula (|45[) 
(dashed line) when .M*© = 5.33 and h=l. The dash-dot-dash horizontal line 
represents the upper limit of the complete sample. 

high values of x 2 are P ar t of expected results. The reduced values of x 2 e d 
conversely give acceptable results. 

From our numerical analysis we draw the following conclusions: 

• The generalized gamma galaxy LF with five parameters, formula (|16[) . 
produces a "better fit" than the Schechter LF in all five bands con- 
sidered, see Table [TJ As an example, in the u* band, \ 2 an d Xred are 
282 and 0.59 against 330 and 0.68 for the Schechter LF. 

• The generalized gamma galaxy LF with four parameters, formula (|2ip . 
produces a "better fit" than the Schechter LF in all five bands con- 
sidered, see Table [2j As an example, in the g* band, x 2 an d X 2 e d are 
747 and 1.263 against 753 and 1.263 for the Schechter LF. Because 
the number of parameters is less, x 2 > AIC and BIC are higher than 
in the case of the generalized gamma galaxy LF with five parameters, 
see Table Q] 

• The LF of galaxies obtained from the generalized gamma PDF for the 
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Fig. 14. The galaxies in the 2dFGRS database with 13.44 < bJmag < 16.94 or 
17950 j^p^ < / < 493844 ; are organized as frequency versus heliocentric 
redshift, (empty stars). The theoretical curves generated by the integral of the 
Schechter function in flux (formula (|3"8")l with parameters as in Table O (full line) 
and by the integral of LF4 as represented by formula (|43p with parameters as in 
column 2dFGRS of Table [5]) (dashed line) are drawn when A^*© = 5.33 and h=l. 
The numerical analysis gives \ 2= 3314 for the Schechter function (full line) and 
X 2 = 2246 for LF4 (dashed line). 



mass of galaxies assuming L oc M d , formula (|28p . produces values of d 
that are slightly greater than one and a " better fit" than the Schechter 
LF in all five bands considered, see Table As an example, in the 
r* band, x 2 and Xred are 2217 and 3.31 against 2260 and 3.36 for the 
Schechter LF. 

The Schechter LF of galaxies can also be improved by incorporating a 
transformation of location, formula (|30p . and the numerical analysis pro- 
duces a " better fit" than the standard Schechter LF in three out of the five 
bands analyzed, see Table HI As an example, in the z* band, \ 2 an d Xred 
are 3133 and 4.25 against 3245 and 4.4 for the Schechter LF. 

One of the four new LFs, more precisely the generalized gamma with 
four parameters, (equation (|18p ). was tested on the 2dFGRS database and 
the analysis of x 2 on the histogram of observed frequencies versus redshift 
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Fig. 15. Range in percent of the observable absolute magnitude as a function of 
the redshift when M max =-1Q, M min — -22, rriL— 19.61 and h=l. 

produces x 2= 6654 for the new LF against x 2 =8078 for the Schechter func- 
tion. 
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